Spin Squeezing via One- Axis Twisting with Coherent Light 
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We propose a new method of spin squeezing of atomic spin, based on the interactions between 
atoms and off-resonant light which are known as paramagnetic Faraday rotation and fictitious mag- 
netic field of light. Since the projection process, squeezed light, or special interactions among the 
atoms are not required in this method, it can be widely applied to many systems. The attainable 
range of the squeezing parameter is £ > S~ 2 ^ 5 , where S is the total spin, which is limited by 
additional fluctuations imposed by coherent light and the spherical nature of the spin distribution. 

PACS numbers: 03.67.Mn, 32.80.-t, 42.50.Lc 



Squeezed spin state (SSS) is one of the non-classical 
states in collective spin system. In SSS, the quantum 
uncertainty of the spins along an axis orthogonal to the 
mean spin vector (AS* 2 ) is suppressed below the standard 
quantum limit (SQL) such as (ASj_) < |(S)|/2, where 
(S) is the mean spin vector, due to an entanglement for- 
mation among the individual spins. The degree of the 
squeezing isjisually evaluated by the squeezing parame- 
ter ( = 2(A5j_)/| (S) |, in terms of the variance to average 
ratioQ. 

For the last several years, SSS has been extensively in- 
terested in not only for precision measurement of a spin 
component but also for the application to the 

quantum infomation[a, |6|, |j| • There have been many pro- 
posals and experiments to realize the spin squeezing of 
atoms. They can be put into three categories as follows 
: (i) Quantum non-demolition (QND) measurement of 
spin via paramagnetic Faraday rotation and spin squeez- 
ing by quantum projection 0, H H H [H El [3 . The 
QND measurement has been already performed by some 
groups for the electronic ground states of atom jol ll llll3| . 
and the squeezed parameter has reached about C ~ 0.7 
for S ~ 4 x 10 7 .11], and C ~ 0.1 for S ~ 1O U 0- Since 
the projection causes the squeezing in this method, the 
degree of the squeezing will be finally determined by the 
performance of the detector, (ii) Quantum-state transfer 
from squeezed light to spin [j| 0, 0, 0, 0, ■ One type 
is based on the complete absorption of squeezed vacuum, 
and has been experimentally demonstrated for the elec- 
tronic exited states of atom (( ~ 0.97 for S ~ 5 x 10 7 ) 
EH- Another type is based on the stimulated Raman 
adiabatic passage 0, 0, E El • Since the squeezed light 
is the source of spin squeezing in these methods, the de- 
gree of the squeezing will be finally determined by the 
quality of the squeezed light, (iii) Special systems to 
induce nonlinear interactions among the individual spins 
such as Bose- Einstein condensates [l8L fl9l Eof . cold atoms 
in optical latticejil]], atoms in optical cavity [22lE^.E4l|. 
They are not easy to prepare and difficult to operate after 
squeezed. 



In this paper, we propose a new method to realize the 
spin squeezing, which can not be put into any of the three 
categories. Our method does not rely on the projection 
by the measurement, use of squeezed light, and the spe- 
cialities of the systems. Instead, the new method only 
requires a coherent light pulse and a few linear optics, so 
it can be widely applied to many systems. It should be 
noted that a recent electronic archive by K. Hammerer et 
aL|25| includes another proposal of an unconditional spin 
squeezing with coherent light. 

Our method is based on the interaction between atoms 
and off-resonant light, whose interaction Hamiltomian 
takes a form 



H = aJ z S z 



(1) 



where a is a real constant, and the z-axis is set paral- 
lel to the wave vector of the light. S is the summation 
over the individual spin, which obeys the usual commu- 
tation relation of angular momenta [S, S] = iS. J is 
quantum- mechanical Stokes vector of light, which also 
obeys the usual commutation relation of angular mo- 
menta [J, J] = iJ. For a light pulse with the dura- 
tion T propagating in free space, J can be written as 
J x = \ Jq \a+a-+aLa + )dt, J y = i ' (a\_a_ — a?_a + )dt, 
J z = | (a^a^- — o^_a )dt, where a± is the annihilation 
operators of a± circular polarization mode, respectively 
[26j . The interaction of Eq. represents the addition of 
the phase difference for a± light, which causes the rota- 
tion of the polarization plane for linear polarization at the 
angular frequency aS z /2, known as paramagnetic Fara- 
day rotation. It also represents the spin rotation around 
the z-axis at the angular frequency aJ z , known as ficti- 
tious magnetic field of light [27] . If we are able to apply a 
light pulse whose J z is proportional to S z as a fitctitious 
magnetic field, the collective spin will nonlinearly rotate 
at angular frequencies proportional to Sz, whose evolu- 
tion will be similar to one-axis twisting This is the 
basic idea of our proposal. 

To design such an interaction, we propose a system 
illustrated in Fig|T] Initially a light pulse is lin- 
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FIG. 1: System of our proposal. A linearly polarized light 
pulse passes through an atomic ensemble and the polariza- 
tion plane is rotated. The rotation angle is proporitonal to 
S z and converted to the circular polarization components af- 
ter passing through the A/8 plate twice. When the pulse 
passes through the atomic ensemble again, the pulse induce a 
nonlinear rotation to the atomic ensemble around the z-axis 
as a fictitious magnetic field. See text for details. 



early polarized along the rc-axis and contains 2J(3> 1) 
photons as an average. Atoms \ips) ar e spin-polarized 
along the £-axis and contains total spin S. The ligh t 
is weakly focused to match the atomic ensemble |26l . 
The averages of the Stokes components is then ( J x ) = J, 
(J y ) = (J z ) = 0, and the averages of the collective spin 
vector is (S x ) = S, (S y ) = (S z ) = 0. Since the light 
pulse is a strong coherent state, we can approximate the 
commutation relation as [J y ,J z ] — iJ |26(. Firstly, a 
light pulse passes through the atoms and the polariza- 
tion plane is then rotated. We call it "the first interac- 
tion" , whose interaction time is labelled as t\ . The Stokes 
vector becomes j( FI ) = e ltlH 3e~ ltlH , whose y compo- 



(FI) 



J-u 



v — u v ~ 

(FI) 



OitiJS z , 

becomes 



nent is approximately written as J, 
for at\S z <C 1. Since the average of the J y 
{^AJy Wj) — a t\JS Zl we can say that the information 

( FI) 

of S z is copied and held on J y ' as a Faraday rotation an- 
gle. We note that S z is conserved because the interaction 
of Eq. satisfies the back-action evasion (B AE) condi- 
tion of [S Z ,H] = 0. Secondly, the pulse passes through 
twice the A/8 wave plate by the totally retroreflecting 
mirror. As a result, A/4 phase difference is induced be- 
tween the two orthogonal modes of linear polarization. 
We call it "the local operation" for the light. The Stokes 
vector becomes j( LO ) = e^l 2 )^ j(Fi) e -i(*72) w h se z 

component is j| LO ' > = J y FV> . We can say that the infor- 
mation of S z is shifted from Jy FI ' to ji LO \ converting 
the angle of the polarization plane to the photon num- 
ber difference of the o± modes. Thus, the required light 
is achieved whose J z is approximately proportional to 
S z . Finally, the pulse passes through the atomic ensem- 
ble again. We call it "the second interaction", whose 
interaction time is labelled as t^- The interaction Hamil- 
tonian of the second interacton can be roughly written 
as i/( SI ) ~ aji h °' ) S z oc S z , which takes a form similar to 
the one-axis twisting Hamiltonian x$ z Q- Thus, we can 
expect that the spin state becomes SSS after the second 
interaction. 



Next, we derive the density operator of the spin af- 
ter the second interaction to calculate the properties 
of the spin state obtained by this method. The ini- 
tial density operator of the whole system can be written 
as p s j = Ps <S> \ipj)(ipj\, where p s = \ips)(ips\- After 
the second interaction, it becomes psj = UpsjU^ where 
U _ e -it 3 H e -Hir/2)J xe -itiH _ Tne re duced density oper- 
ator representing the spin state after the second interac- 
tion ps can be written as ps = Trj(psj), where Trj is the 
partial trace for the light. For convenience, we consider 
the set of eigenstates for S 2 and S z , say \S,M), where 
S 2 \S,M) = S(S + 1)\S,M) and S Z \S,M) = M\S,M). 
The matrix elements take a form 

(S, M\p s \S, M') = <j mm 1 {S, M\ PS \S, M'), (2) 

°MM* = e -^^-M') 2 /2 e -MM*~M")/2^ (3) 

where we have set p, = {ocb\){<xt2)J and p' = ((ati) 2 + 
{ott2) 2 )J '/2. If t\ = ti then p = pi . Since the atoms \ips) 
are polarized along the x-axis, the matrix elements of ps 
can be written as 

1 / 25 \ 1/2 / 25 ^ 1/2 



(S,M\ PS \S,M') = ^ S - M ) ^; M/ j . (4) 

In the following discussions, we use the experessions of 
Eq.© and Eq.|@}. Wc note that the ideal one-axis 
twisted state [l| corresponds to the case of pi = 0. 

To know how uncertainties evolve, we calculate 
the quasiprobability distributions (QPD), which is de- 
fined as Q{6,4>) = (6,<j)\ps\0,4>), where \6,<f>) = 



i6S v 



S, S) is a spin state polarized along the di- 
rection whose polar and azimuth angles are 9 and <fi, 
respectively 0. The results of the calculations in the 
case of S = 20 are shown in FigEl for the initial spin 
state (a), the spin state after the first interaction (b), 
and that after the second (c). The initial spin state is 




FIG. 2: State evolutions expressed as the quasiprobability 
distribution for S = 20. The value of QPD for (9, (j)) direction 
is represented by the gray scale on the unit sphere, which is 
normalized by the maximum value, (a) The initial spin state, 
(b) The spin state after the first interaction, where we have 
set (ati) 2 J/2 — 0.1 and = 0, in other words, p = and 
p! = 0.1. (c) The spin state after the second interaction, 
where we have set (crfi) 2 J/2 = (at 2 ) 2 J/2 = 0.1. in other 
words, p, = p! — 0.2. The spin squeezing is realized along the 
z'-axis. 

isotropically distributed along the x-direction as is shown 
in Figl^a). After the first interaction, the distribution 
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is a little broadened along the y direction, as Fig^b) 
indicates. This is explained by additional fluctuation 
imposed by coherent light. In fact, the y components 
after the first interaction is approximately written as 
s m = e i tl H Sye -iHH ~ s y + at 1 J z S x for ah J z < 1. 
Since S z is the BAE variable, the distribution along the 
z direction does not change at all. After the second in- 
teraction, the distribution looks twisted around the z- 
axis and squeezed along the z'-axis, as Fig^fc) indi- 
cates. Although not clear from the figure, the distri- 
bution is also broadened along the y-axis as in the case 
after the first interaction. In fact, the y component af- 
ter the second interaction is roughly written as S y ~ 
e it 2 H^ s m e -it 2 H^ „ Sy +f,S z S x + (atiJ z +od 2 J y )S x . 
By these additional fluctuations imposed by coherent 
light, the spin state after the second interaction is dif- 
ferent from the ideal one-axis twisted state 0. The ad- 
ditional fluctuations would be reduced by use of a polar- 
ization squeezed light pulse whose squeezed component 
is t\J z -\-tiJ y , approaching the ideal one- axis twisting in- 
teraction of fi' — > Q. We mention that the additional fluc- 
tuations by light in the method of Fief. [2^1 are imposed 
both on the y and z components, while the z component 
is squeezed. Therefore, the squeezing parameter does not 
become small in that scheme. 

From Eq.@ and Eq.Q, we can derive the averages 
and variances of the spin components. The averages 
can be calculated as (5a,) = 5e~ M / 2 cos 2S ~ 1 (fi/2) and 
(S y ) — (S z ) — 0, where S represents the spin operator 
after the second interaction. They indicates that the ori- 
entation of the mean spin vector remains (w/2,0) direc- 
tion or the x-axis, as is shown in Fig^c). To character- 
ize a elliptical distribution around the x-axis, we define 
the minor and major axis, say z' and y' , respectively, 
as is shown in FigE^c), so that the variances of whose 
components (AS 2 ,) and (AS 2 ,) give the minimum and 
maximum on the y — z plane, respectively. The variances 
can be calculated as (A5 2 ) = S 2 - (S x ) 2 - S{S -1/2) A/2 
and 



(ASl 



S 
2 



SS_ 
2~ 



1/2 



(A ± y/A 2 + B 2 



(5) 



where we have set A = 1 



-V 



cos b fi, and B 



46-^/2 gin^/2) cos 2S " 2 (/i/2). We can also calculate 5, 
which is an angle between the directions of the z'- and 
z-axes, or the y'- and y-axes, as is shown in Fig|2Jc), 
and obtain 5 = arctan(i3/^4)/2. For S > 1 and 5 _1 < 
/i ~ // <C we find the approximate value of the 

variance of the z' component 



(ASl)c^ 7 



2 V7 2 +7' 3 



1 + 



(6) 



where we have set 7 = Sfi/2, 7' = Sfi'/2 and j3 = Sfi 2 /4. 
Also we find (S x ) ~S(l-/3). 



To examine the dependence on the interaction strength 
ati,at2 and the input photon number 2 J, we plot the 
variances of the y' components and the z 1 components as 
a function of fi(— fi') in Fig|3fa). We also plot the ap- 
proximate value for the z' components written as Eq.©. 
It is clearly known that the variance of the z' component 
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FIG. 3: (a)Variances of the z' (circle below 1) and the y' 
(circle above 1) components for S = 20 as a function of 
M (Mhaif = 0.117, (j, m i n — 0.236). They are normalized as 
2(A£y)/S and 2{AS y ,)/S, respectively. We also plot the ap- 
proximate value of 2{AS^,)/S (solid line), (b) Value of /ihaif 
(triangle) and fi m in (square) as a function of S, which are the 
required values of fi to obtain the half variance of the SQL 
and the attainable minimum variance, respectively. We also 
plot the approximate solutions of /ihaif (dashed line) and fi m ixx 
(solid line). To obtain both (a) and (b), we have assumed 
M = M- 

is reduced for small fi, minimized at an optimal value 
of fi, and becomes large for large fi. It means that too 
strong interaction or too large photon number deterio- 
rates the squeezing. This is explained by the spherical 
nature of the spin distribution, and in fact, the variance 
of the y' component is almost saturated at the largest 
value of 5 2 /2 for large fi, which was entirely ignored in 
the analysis in Ref. 25] . As a typical value of fi, we in- 
troduce /ihaif as the value of ft to attain (AS 2 ,) = 5/4, 
the half variance of the SQL. We also introduce /i m i n as 
the value to attain the minimum of the (AS 2 ,). We plot 
the numerical solutions of /ihaif and /i m in in Fig[3| (b) for 
the case of fi = ft! . One can see that both /ihaif and 
A*min become small as 5 increases but they obey differ- 
ent power laws. From Eq.@), we find /ihaif — 25 _1 and 
/i m ; n ~ 2(3/2) 1 / 5 5 -3 / 5 . We show these approximate so- 
lutions in Fig|3|(b), which are in good agreement with the 
numerical ones. We also find that the squeezing parame- 
ter at /i = fi 1 = /i min becomes Cmin — (2/3) 1 / 5 5 -2 / 5 . We 
note that it is slightly worse than (l/3) 1 / 3 5 -2 / 3 , which is 
the squeezing parameter for the ideal one-axis twisting, 
due to the additional fluctuations imposed by coherent 
light, as is mentioned above. 

Finally, we discuss the feasibility of our method. In 
the following, we consider the case that the shape of the 
light pulse is a square wave with its peak power P and 
pulse duration T. As in Ref. [3], we assume A ^> T, 
il <C A, and rT <§C (5/i) _1 <C 1, where A represents the 
detuning from the resonance frequency, T the full natural 
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linewidth at half maximum of the transition, fl the Rabi 
frequency, and r the photon scattering rate |2s|. After 
some calculations, we obtain fi = fi' = rTa / (2ttw 2 ), 
where w represents the beam waist and cxo the photon- 
absorption cross section of an atom, which can be writ- 
ten as o"o = 3Aq/(27t) with the resonance wavelength 
Ao- We note that fi, r and /j, are exactly the same as 
2g^2N p /{cT), 4e a /T, and n 2 /N a in Ref.0, respec- 
tively. The condition to obtain fi > jt%aif, or £ < 1/2, can 
be rewritten as do ^ 8(rT) _1 , where do = 2Sao/(nw 2 ) 
is the optical depth. We also note that this condition is 
the same as k > y/2, similar to that of QND measure- 
ment |26j. Such a condition has been satisfied in several 
systems, such as atoms in a cell laser cooled and 
trapped atom, and so on. The feasibility of our scheme 
also comes from the simple experimental setup depicted 
in Fig^ which is also the great advantage over another 
scheme in Ref.|25j. This suggests efficient squeezing can 
be realized by the current technologies. 

As one ideal example, we consider a ytterbium atoms 
( m Yb) in optical trapji^, lief , which contains S = 
4 x 10 6 . The atom collision and the precession due to the 
stray magnetic field, which causes the transverse relax- 
ation, are well surpressed because it is ultracold fermion 
and has only nuclear spin 1/2 whose gyromagnetic ratio 
is about three orders smaller than paramagnetic atoms 
like alkali metal. From the parameters given in Ref . [30| . 
w = 3fim 30], A = 399nm and T = 2tt x 29MHz, the 
light pulse of /i = 5.4 x 10~ 6 , for example, is obtained by 
setting A = 2tt x 24GHz, P = 17nW, T = 0.24ms, which 
satisfies the assumptions A ^> T, Q = 2n x 21MHz <C A, 
rT = 4.0 x 10~ 3 < (S>) _1 = 4.6 x 10~ 2 < 1. and 
J = 4.0 x 10 6 3> 1. In this case the squeezing parameter 
becomes £ = 0.08. We note that the decay constant of the 
atom number in optical trap is about 4s [3jJ , which is so 
longer than the the pulse duration T that we can treat 
the total spin S as a constant. We also note that the 
length of the atom distribution shuould be adjusted to 
L ~ 70/im to satisfy the condition of irw 2 / (XqL) ~ l|26|. 
which is easy for atoms in optical trap of crossed config- 
uration. 

To avoid the interference between the three steps of 
the first interaction, the local operation and the second 
interaction within one pulse, we can use a pluse train, 
each duration of which is so short that the three steps 
are separable and the repetition rate is so slow that the 
pulse number travelling in the path is at most one. Since 
the each matrix element (S, M\ps\S, M'} would evolve 
like a geometric progression whose common ratio is a mm 1 
for the every pulse as Eq.J2J indicates, and a mm 1 is the 
power of the mean photon number of the each pulse as 
Eq.© indicates, we can say that the same SSS would 
be obtained as long as the total mean photon number 
passed through the atomic ensemble are equal. 
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